Candidate-Based Approaches to Identify Genetic Variation Influencing Type 2 Diabetes and Quantitative Traits

Type 2 diabetes (T2D) is a metabolic disorder characterized by insulin resistance and impaired insulin secretion that affects more than 20 million Americans, although the genetic component of the disorder is largely unknown. Individual genetic susceptibility to type 2 diabetes and other complex traits is the result of variation that is both common in human populations and rare, de novo and inherited mutations. We adopted a diverse set of genetics, genomics and informatics approaches to prioritize candidate genomic regions and variants and perform in-depth, targeted analysis of their contributions to type 2 diabetes susceptibility and related trait variability. Our initial efforts focused on the selection of candidate genes relevant to a complex trait by developing a metric to weight the relevance of functional gene annotations to the known biology of a trait. We used this method to select candidate genes for type 2 diabetes and performed a T2D case-control and quantitative trait association study in 2,335 Finnish individuals from the FUSION study. After follow-up in additional samples, we identified several variants that might contribute to T2D susceptibility. Genomic regions associated with plasma levels of HDL cholesterol and triglycerides were re-sequenced in individuals with trait-extreme values. Our analysis revealed a denser set of common and rare functional target variants including several non-synonymous, 3' UTR, and non-coding SNPs and indels. Finally, we utilized two approaches to identify candidate functional non-coding variants that may directly contribute to trait susceptibility. First, we used Formaldehyde-assisted isolation of regulatory elements (FAIRE) coupled with high-throughput sequencing to identify nucleosome-depleted regions in pancreatic islets. We used islet FAIRE-seq data to identify SNPs associated with T2D that potentially alter islet transcriptional regulation. A SNP in TCF7L2, rs7903146, was located in a FAIRE-seq site and demonstrated allelic differences in islet chromatin openness and enhancer activity, suggesting that it may contribute functionally to T2D susceptibility. Second, we used transcription factor binding site motifs to computationally predict variants that have allelic differences in regulatory activity. Taken together, these results suggest that identifying candidate genomic regions can successfully enrich for variation important for type 2 diabetes and other complex traits

Lipid-Induced Epigenomic Changes in Human Macrophages Identify a Coronary Artery Disease-Associated Variant that Regulates \u3cem\u3ePPAP2B\u3c/em\u3e Expression through Altered C/EBP-Beta Binding

Genome-wide association studies (GWAS) have identified over 40 loci that affect risk of coronary artery disease (CAD) and the causal mechanisms at the majority of loci are unknown. Recent studies have suggested that many causal GWAS variants influence disease through altered transcriptional regulation in disease-relevant cell types. We explored changes in transcriptional regulation during a key pathophysiological event in CAD, the environmental lipid-induced transformation of macrophages to lipid-laden foam cells. We used a combination of open chromatin mapping with formaldehyde-assisted isolation of regulatory elements (FAIRE-seq) and enhancer and transcription factor mapping using chromatin immuno-precipitation (ChIP-seq) in primary human macrophages before and after exposure to atherogenic oxidized low-density lipoprotein (oxLDL), with resultant foam cell formation. OxLDL-induced foam cell formation was associated with changes in a subset of open chromatin and active enhancer sites that strongly correlated with expression changes of nearby genes. OxLDL-regulated enhancers were enriched for several transcription factors including C/EBP-beta, which has no previously documented role in foam cell formation. OxLDL exposure up-regulated C/EBP-beta expression and increased genomic binding events, most prominently around genes involved in inflammatory response pathways. Variants at CAD-associated loci were significantly and specifically enriched in the subset of chromatin sites altered by oxLDL exposure, including rs72664324 in an oxLDL-induced enhancer at the PPAP2B locus. OxLDL increased C/EBP beta binding to this site and C/EBP beta binding and enhancer activity were stronger with the protective A allele of rs72664324. In addition, expression of the PPAP2B protein product LPP3 was present in foam cells in human atherosclerotic plaques and oxLDL exposure up-regulated LPP3 in macrophages resulting in increased degradation of pro-inflammatory mediators. Our results demonstrate a genetic mechanism contributing to CAD risk at the PPAP2B locus and highlight the value of studying epigenetic changes in disease processes involving pathogenic environmental stimuli

Common DNA sequence variation influences 3-dimensional conformation of the human genome.

BACKGROUND:The 3-dimensional (3D) conformation of chromatin inside the nucleus is integral to a variety of nuclear processes including transcriptional regulation, DNA replication, and DNA damage repair. Aberrations in 3D chromatin conformation have been implicated in developmental abnormalities and cancer. Despite the importance of 3D chromatin conformation to cellular function and human health, little is known about how 3D chromatin conformation varies in the human population, or whether DNA sequence variation between individuals influences 3D chromatin conformation. RESULTS:To address these questions, we perform Hi-C on lymphoblastoid cell lines from 20 individuals. We identify thousands of regions across the genome where 3D chromatin conformation varies between individuals and find that this variation is often accompanied by variation in gene expression, histone modifications, and transcription factor binding. Moreover, we find that DNA sequence variation influences several features of 3D chromatin conformation including loop strength, contact insulation, contact directionality, and density of local cis contacts. We map hundreds of quantitative trait loci associated with 3D chromatin features and find evidence that some of these same variants are associated at modest levels with other molecular phenotypes as well as complex disease risk. CONCLUSION:Our results demonstrate that common DNA sequence variants can influence 3D chromatin conformation, pointing to a more pervasive role for 3D chromatin conformation in human phenotypic variation than previously recognized

Identification of a Regulatory Variant That Binds FOXA1 and FOXA2 at the CDC123/CAMK1D Type 2 Diabetes GWAS Locus

Many of the type 2 diabetes loci identified through genome-wide association studies localize to non-protein-coding intronic and intergenic regions and likely contain variants that regulate gene transcription. The CDC123/CAMK1D type 2 diabetes association signal on chromosome 10 spans an intergenic region between CDC123 and CAMK1D and also overlaps the CDC123 3′UTR. To gain insight into the molecular mechanisms underlying the association signal, we used open chromatin, histone modifications and transcription factor ChIP-seq data sets from type 2 diabetes-relevant cell types to identify SNPs overlapping predicted regulatory regions. Two regions containing type 2 diabetes-associated variants were tested for enhancer activity using luciferase reporter assays. One SNP, rs11257655, displayed allelic differences in transcriptional enhancer activity in 832/13 and MIN6 insulinoma cells as well as in human HepG2 hepatocellular carcinoma cells. The rs11257655 risk allele T showed greater transcriptional activity than the non-risk allele C in all cell types tested. Using electromobility shift and supershift assays we demonstrated that the rs11257655 risk allele showed allele-specific binding to FOXA1 and FOXA2. We validated FOXA1 and FOXA2 enrichment at the rs11257655 risk allele using allele-specific ChIP in human islets. These results suggest that rs11257655 affects transcriptional activity through altered binding of a protein complex that includes FOXA1 and FOXA2, providing a potential molecular mechanism at this GWAS locus

The genetic architecture of type 2 diabetes

The genetic architecture of common traits, including the number, frequency, and effect sizes of inherited variants that contribute to individual risk, has been long debated. Genome-wide association studies have identified scores of common variants associated with type 2 diabetes, but in aggregate, these explain only a fraction of heritability. To test the hypothesis that lower-frequency variants explain much of the remainder, the GoT2D and T2D-GENES consortia performed whole genome sequencing in 2,657 Europeans with and without diabetes, and exome sequencing in a total of 12,940 subjects from five ancestral groups. To increase statistical power, we expanded sample size via genotyping and imputation in a further 111,548 subjects. Variants associated with type 2 diabetes after sequencing were overwhelmingly common and most fell within regions previously identified by genome-wide association studies. Comprehensive enumeration of sequence variation is necessary to identify functional alleles that provide important clues to disease pathophysiology, but large-scale sequencing does not support a major role for lower-frequency variants in predisposition to type 2 diabetes

Single-trait and multi-trait genome-wide association analyses identify novel loci for blood pressure in African-ancestry populations.

Hypertension is a leading cause of global disease, mortality, and disability. While individuals of African descent suffer a disproportionate burden of hypertension and its complications, they have been underrepresented in genetic studies. To identify novel susceptibility loci for blood pressure and hypertension in people of African ancestry, we performed both single and multiple-trait genome-wide association analyses. We analyzed 21 genome-wide association studies comprised of 31,968 individuals of African ancestry, and validated our results with additional 54,395 individuals from multi-ethnic studies. These analyses identified nine loci with eleven independent variants which reached genome-wide significance (P \u3c 1.25×10−8) for either systolic and diastolic blood pressure, hypertension, or for combined traits. Single-trait analyses identified two loci (TARID/TCF21 and LLPH/TMBIM4) and multiple-trait analyses identified one novel locus (FRMD3) for blood pressure. At these three loci, as well as at GRP20/CDH17, associated variants had alleles common only in African-ancestry populations. Functional annotation showed enrichment for genes expressed in immune and kidney cells, as well as in heart and vascular cells/tissues. Experiments driven by these findings and using angiotensin-II induced hypertension in mice showed altered kidney mRNA expression of six genes, suggesting their potential role in hypertension. Our study provides new evidence for genes related to hypertension susceptibility, and the need to study African-ancestry populations in order to identify biologic factors contributing to hypertension

Allele-specific NKX2-5 binding underlies multiple genetic associations with human electrocardiographic traits.

The cardiac transcription factor (TF) gene NKX2-5 has been associated with electrocardiographic (EKG) traits through genome-wide association studies (GWASs), but the extent to which differential binding of NKX2-5 at common regulatory variants contributes to these traits has not yet been studied. We analyzed transcriptomic and epigenomic data from induced pluripotent stem cell-derived cardiomyocytes from seven related individuals, and identified ~2,000 single-nucleotide variants associated with allele-specific effects (ASE-SNVs) on NKX2-5 binding. NKX2-5 ASE-SNVs were enriched for altered TF motifs, for heart-specific expression quantitative trait loci and for EKG GWAS signals. Using fine-mapping combined with epigenomic data from induced pluripotent stem cell-derived cardiomyocytes, we prioritized candidate causal variants for EKG traits, many of which were NKX2-5 ASE-SNVs. Experimentally characterizing two NKX2-5 ASE-SNVs (rs3807989 and rs590041) showed that they modulate the expression of target genes via differential protein binding in cardiac cells, indicating that they are functional variants underlying EKG GWAS signals. Our results show that differential NKX2-5 binding at numerous regulatory variants across the genome contributes to EKG phenotypes